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A MATHEMATICAL FORMULATION OF THE SCOLE CONTROL PROBLEM: PART 1 


SUMMARY 


This is the first part of a multi-part report on a mathematical formula- 
tion of the SCOLE control problem, and deals primarily with various techniques 
of solution of the partial differential equations with delta-functions on 
the boundary, the most comprehensive of which is the formulation as a non- 
linear abstract wave equation which clarifies the concept of modes and leads 
immediately to a linear feedback law for stability augmentation. We also 
obtain a "closed-form solution" in a special case. 


1 . INTRODUCTION 

This report is in several parts. In this the first part we deal mainly 
with techniques of solution of the partial differential equations with delta- 
functions on the boundary, as formulated in [1] . Several methods of solution 
are described. The most comprehensive of these is the formulation as a non- 
linear abstract wave-equation, which clarifies the concept of modes and allows 
us to develop a linear feedback that can stabilize the system, in the sense 
that the total energy cannot increase. We are also able to obtain a "closed- 
form solution" in a special case, using a "boundary input" approach following 
[3]. Finally, the mathematical formulation of the problem should help in the 
digital computer simulation of various control laws. 


A/S5~XZ c iO^ 
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While this report closely follows [1] in notation, the beginning section, 
Section 2, restates the Beam Equations in such a way as to clarify the delta- 
function formulation in [1]. Section 3 outlines a Green f s Function approach 
to the solution. Section 4 shows how the 6-functions can be removed, and 
replaced by time-varying boundary conditions. Section 5 shows how the Laplace 
Transform technique could be used and shows its equivalence to the Green’s 
Function method. Finally Section 6 outlines the "Boundary Input" method patterened 
on the development in [3] . After specializing the boundary inputs to the spe- 
cified boundary dynamics (in terms of ordinary differential equations) we go on 
to formulate an "abstract wave equation" in Section 8. Using this formulation, 
we are able to show that it is possible to stabilize the system (in terms of 
decreasing the total energy of the system) by a linear feedback. Moreover 
we also indicate how to calculate the modes, as well as a modal expansion of 
the solution leading to a nonlinear integral equation to take care of the 
nonlinearity introduced in the rotational motion. Finally in Section 9 we 
indicate one explicit solution using the "boundary input" technique of Section 6. 
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2. THE BEAM BENDING EQUATIONS 


U(t, s) = 


u «{, (t ’ s) 

u 0 (t,s) 

u ^ ( t, s ) 


0 < t, -°° < s < 


U(t,s) = 


s < 0, s > L . 


Superdots will denote derivatives with respect to time t, 0 < t. Primes 
will denote derivatives with respect to the space variable s. The basic 
beam deflection equations are: 


Roll Beam Bending Equation 


PA iijCt.s) + EI^J"(t,s) 


I S(s - s n> + I «•(.-.„) , 


0 < t, -CO < s < 


0 < s < L . 
n 


( 2 . 1 ) 


The delta -function and the delta-function derivative are to be interpreted as 
generalized functions. Hence for any C*" function h(s), < s < we 


L 4 4 

/ <PAu.(t, S ) + EI.u""(t,s)> h(s) ds = S „(t)h(s) + I g, „(t)h'(s > . 


( 2 . 2 ) 


Pitch Beam Bending Equation 


•u e (t,s) + EI e u""(t,s) - I f 6>ti (t) 6(s-s n ) + l g 0>n (t) «’(s-s n ) 


0 < t, -CO < S < 


0 < s < L 
n 


(2.3) 
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The interpretation of the delta-functions is the same as in the case of the 
previous Roll equations, so need not be repeated. 

Yaw Beam Torsion Equation 

4 

= I ^( s-s n ) (2.4) 

0<t, -^<s<°°, 0 < s < L 

n 

Once again, the interpretation of the delta -function is the same as before. 


3. METHOD OF SOLUTION I: GREEN'S FUNCTIONS 

Roll Beam Bending Equations 

We begin by constructing the Green’s function for the homogeneous equations: 
PA u^(t,s) + EI^uJ"(t,s) = 0, 0 < t, 0 < s < L 

with the boundary conditions ("free-free") 
u^(t,04-) = u"(t,L-) = 0 
uj"(t,0+) = u”’(t,L-) = 0 

As is well-known [4], this can be done with the aid of the eigenfunctions 
(*) which satisfy 

“4 V s) = k n V s) ’ 0 < s < L 

as 

and (to satisfy the boundary conditions) is of the form: 

♦n* 8 ) = A n (cosh k n s + cos k s) + B (sinh k s + sin k s) 

ii il n n n n 
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where A , B are such that 
ir n 


and, of course, 


The k^ satisfy 


/ <f> (s) 2 ds - 1 


/ t n w 4> m (s> ds - 


5 m 

n 


cosh k n L cos k^L = 1 


It is well-known that [4] 


k„ = 0 


Letting 


k ^ 
n 


Y, 


(n+l)ir 

2L 


ft 

PA * 


the Green 1 s function is defined by 


00 sin y,k t 

G(j) ( t » s , s ' ) = l n 


0 


Y . k 

<P n 


2 — 4> n (s) <J> n (s') , 0 < t, OS s,s' < L . 

(3.1) 


We may now apply Duhamel’s principle [5] to obtain the solution to the 
non-homogeneous Roll Beam Bending Equations (2.1). Thus we have 
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u 4) (t >S ) 


t 4 L 

I l I G,(t-o; s; S')(f. (a) 6 (s' -s') + g. (o)S' (s' -s')) da ds’ 


0 10 


n B 4>»n 


+ y a sin (y,k 2 t + 3 ) 4> (s) 
k n '<p n n' r n ' 


/ I yt-o; =; s n ) f^to) do + ! I G’<t-o; S, s n ) g^Co) d0 


+ V a sin (y,k 2 t + 3 ) $ (s) 
k n n n' Y n w 


(3.2) 


where 


G ^ t » s » s n ^ 


gjT- G (t,s,s*) 

s =s 


The term outside the integral in (3,2) takes care of the initial condi- 


tions: 


u <f>(o;s)» u^(o,s> 


for suitable choice of { a and C^}, as is well known [4]. For simplicity, 
we shall often set the initial conditions to be zero. 


Pitch Beam Bending Equations 


In a similar fashion, defining 


r e pa 


and the Green’s function 


Gq (t , s , s ’ ) - 


00 sin v A k t 

1 — v- 
0 


2 ~ <f> n (s) <() n (s’) 


with ^ as before » can express the solution to the Pitch Beam Bending 
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Equations (2.3) as: 
t 4 


t 4 


u 0 (t,s) = f l G (t -o; s, s ) f B (a) do + J £ G'(t-a; s, s ) g fl (a) da 

Q]_ u 01 no,n 


+ y b sin (y Q k^t + 0 ) cf> (s) 
L n 0 n n n 


0 


(3.3) 


Yaw Beam Torsion Equations 

We begin by constructing the Green's function for the homogeneous equations 


p yy t,s) ■ G v? (t,s) = 0 


0<t, 0 < s < L 


with the boundary conditions: 


u ' ( t , 0+) = u'(t,L-) = 0 


(and zero initial conditions). It is well known [4] that this can be done 
with the aid of the eigen functions 4i^(*) which satisfy 


^ (s) = 


ds 


2 r n 


-X ip (s) 
n n 


n 


Let 


Y 




nu 

L ’ 


ft 

PI * 


n = 0, ±1, ±2 


Then the Green's function G^(t,s,s') is given by 


i 00 sin A Y.t 

1 v n'lp , , 

- y L — ~ — cos A s cos A s 
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By Duhamel's principle [5], as before, we can then express the solution to 
the non-homogeneous. Yaw Beam Torsion Equation (2.4) as: 


t L N 

= / / Gw/t-o; s; s’) [g, (a)«(s'-s) dads' 

o o y Vl n j 


^(t.s) 


4 t 


= l I yt-o; s > O g,„ „<<*> <1° 


1 0 


n 


and the general solution allowing for nonzero initial conditions is then 


4 t 


yt,s) 


I / G ^-o; s ‘> s n ) g^ >n (a) da + ^ c n sin (X^t + cos A n s 

(3.4) 


4. METHOD OF SOLUTION II: TIMEVARYING BOUNDARY CONDITIONS 

In this section we shall reformulate the basic equations (2.1) - (2.4) removing 
the delta-functions and replacing them by time-varying boundary conditions. 

Roll Beam Bending Equations 

We begin with the Roll Beam Bending equations (2.1), in the form (2.2). 

Choose h(* ) such that 

h(s) =0 s > e. 

Remembering that 

s^ = 0 

we obtain 

£ 

/ (PAyt,s) + EIjU^'U.s)) h (s) ds = 


ff>,i (t) h(0) + g^yt) h*(o) . 
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We integrate the second integral on the left by parts, twice in succession: 


/ u""(t,s) h(s) ds = u"*(t,e) h(e) - u"’(t,0-) h(0) - u"(t,e)h*(e) 


+ uj(t,0-) h* (0) + / uj(t,s) h"(s) ds . 


Next we let e go to zero. Then we obtain in the limit: 


E y u <f" (t ’ 0+) - U J’ (t. 0 -)) h <°> + El^CuJCt.O-) - uj(t,0+)) h»(0) 


= f <j),l (t) h(0) + 8^ 3.(0 h ' (0) 


Since h(*) can be chosen arbitrarily we obtain that: 


EI 4 ,(u“ (t ,0 + ) - - f M (t) 


E y„i' <t - 0+> " 


Since, by definition, 


uj(t,0-) = uj*(t,0-) = 0 


we obtain 


E W (t > 0+) ' 




In an entirely similar manner, noting that u^" (t,L+) = 
obtain: 


uJ(t,L+) 


= 0, we 


E Vi ,,(t ’ L - ) ■ ‘V (t) 


E Vf (t ’ L - ) - • J 
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Also 

VS" (t ’ S 2 + > ■ f K2 (t) ' 

t 

VS" (t<a 3 +) - u i ,(t,s 3 - )) ' S, 3 (t) • - 

Hence we may replace (2.1) by: 

PA ^(t.s) + EI^u^"(t,s) = 0, 0 < s < L, 0 < t 

with the time-varying boundary conditions: 

■Vi" ■ -S,4 (t) 

VS" (t>s 2 +) - S"^"” ■ S,2 (t) 

VS'^’V’ -S'< t> 0 3 - )) - -S, 3 (t) 

E Vi' (t ’ 0+) ■ - s *,i (t) 

E V» (t - L - > " e 4,,4 (t) 

and any initial conditions in addition that may be imposed: 

u ( j ) (0+, s ) given 

u<j)(0+,s) given 


Pitch Beam Bending Equations 

In an entirely similar manner, we 

PA u Q (t,s) + EIqUq' "( t, s) 


may replace 


(2.3) by: 


(4.3) 


(4.4) 


(4.4a) 


0 


0 < s < L, 


0 < t 


(4.5) 
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with the time-varying boundary conditions 


ElgUg* (t,0+) 




EI e u£' (t,L-) 


- £ 0 , 4 (t) 


ei q( u q' (t,s 2 +) -u^"(t,s 2 -)) 


£ e,2< £ > 


EI 0 (u£' (t,s 3 +) -u^"(t,s 3 -)) 


f 0,3 (t) 


EI q Uq ( t,0+) = 


8 0,l (t) 


EI 0 Ug(t,L-) = 


S 6,4^ 


and any initial conditions imposed 


u fl (0 +, s) given 


u fl (0+, s) given 


Yaw Beam Torsion Equations 


Analogously we may replace (2.4) by: 


PI ^(t, s ) - GljU^'Ct.s) = 0, 0 < s < L, 0 < t 


with the time-varying boundary conditions 


G Vi (t ’ 0+> = - g 4»,i (t) 


'Vi"' 1 -’ ■ + %i (t> 


and any initial conditions imposed: 


u,(0+, s) given 


y°+, s ) 


given 
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5. METHOD OF SOLUTION III: LAPLACE TRANSFORMS 

In this section we shall develop the Laplace Transform technique of 
solution of the Beam equations using the boundary-condition formulation of 
the previous section. We shall take the initial conditions to be zero. Let 


— » — 2 £ 

u(z,s) = J e u(t,s) dt , Re z > 0 
0 


Roll Beam Bending Equations 

The Roll Beam equations (4.4) may be Laplace transformed to yield: 

z 2 U^z.s) + Y 2 u£"(z,s) = 0 , 0 < s < L . 


As is well known, the general solution to this equation for each fixed z 
can be written 


u^(z,s) 


A^(z) cosh 


jr 


s + A 2 (z) cos-J-2- s + A„(z) sinh \— 


+ A^(z) 



(5.1) 


The functions A^(z), A 2 (z) , A^(z) and A^(z) are then to be determined 

by the imposed boundary conditions (4.4a). Thus let 


<f>,n 


*<J>»n 


(z) 

* — .0 

tl 

e Zt f , (t) 

(J),n 

dt 


00 



(z) 

= I 

0 

-zt . s 

6 H.n (t> 

dt 


First, for simplicity, we omit the "proof -mass" conditions at s 2 and s.j. 


Then we have 
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V (z) = E Vi ,,(z ’ 0+) 


3 



+ A„(z) sinh L - A, (z) sintf^- L) 

° ,T <}) h f »<j> j 

There are four linear equations, which can be solved to determine the unknown 
functions. A^z), A 2 (z), A^(z), A^(z), in terms of f^ ^z), f^ ^(z), 

^( z ) anc ^ 4 ( 2 :) • The solution will be unique, as can easily be seen 
(by the non-singularity at the determinant) . 

To include the proof-mass conditions at s 2 and s^ we must subdivide the 
interval [0,L] into three: [0,s 2 ], [s^s^] and [s 3 ,L]. In each of these 

intervals we will express the solution in terms of the coefficient functions 
which will differ now from interval to interval. We omit the details. 

Pitch Beam Bending Equations 

Since the Pitch Beam Bending Equations (4.5) so completely parallel the Roll 
Beam Bending Equation, we shall not need to go into the details of Laplace 
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Transformation. Thus, if we omit the "proof-mass" equations, we can write: 


u 0 (z,s) = B 1 (z) cosh^^- s + B 2 (z) cos-^-s + 


+ B, (z) sin-J— s 

H ,y e 


B,(z) sinh v-p- s 

,y e 


(5.2) 


where the functions B-^z), B 2 (z), B 3 (z), B^(z) are determined as linear 

transformations of: 




* / e ~ Zt £ n JO dt 


-zt 


8 0,n (z) = / e “ „(t) d t 


3 0,n 


just as in the Roll Beam Bending case. To include "proof-mass" conditions, we 
proceed also exactly as in that case. 

Yaw Beam Torsion Equations 

Here Laplace transforming (4*6), we have: 


z\(z,s) 


VJ (z > s) = 0 ’ 


0 < s < L, Re Z > 0 


The general solution is of course 


u^(z,s) = (^(z) cosh ^ + C 2 (z) sinh ^ . 

The substitution of the boundary conditions yield: 


-^,i <z) = G Vt <2 ’ 0+) ' GI 4 ^] c 2 <z) 


S *,* 00 ■ GI . 




— | C. (z) sinh — + C 0 (z) cosh - L 


Y 
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where 


^,l (z) = / e_Zt ^,l (t) dt 

oo 

= / e_Zt g^,4 (t) dt • 


The boundary conditions are simple enough so that we can write the explicit 
solution: 


u^(z,s) = 


^,1 


(z) cosh — (L-s) 


Y 




. , zL 
z sinh — 


Y, 




cosh f 

. , zL 
z smh — 


(5.3) 


where of course we have set all the initial conditions to zero. It is not 

difficult to see that (5.3) is exactly the Laplace transform of the time- 

domain solution expressed by (3.2). We omit the details, except to remark that 

zL 

we need to use the fact that sinh — is an entire function with zeros at 


z 






= X 


Y 


n 




We can make the same statement for the other Beam equation solutions: that 
(5.1) is precisely the Laplace transform of the time-domain solution (3.2) 
and similarly also (5.2) is the Laplace transform version of (3.3). 


6. METHOD OF SOLUTION IV: BOUNDARY INPUTS 

In this technique of solution for the Beam equations, we seek to express 
the solution as the sum of twD . We begin by considering the case where the 
"proof-mass" conditions are omitted. 
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Roll Beam Bending Equation 

Let yt,s) denote the solution of the non-homogeneous equations 
PA ii°(t,s) + EI^u°""(t,s) = z(t,s) , 0 < s < L, 0<t 

with the zero "free— free” end conditions: 


u°"’(t, 0 +) = u°"'(t,L-) = 0 


u°"(t,0+) = u°"(t,L-) = 0 . 


Let 

D^(t,s) = a 3 (t)s 3 + a 2 (t)s 2 + a 1 (t)s + a Q (t) . 

Note that 


d 

— 4 D, (t,s) = 0 , 0 < s < L 

9s y 


where the functions a ± (t), i = 0 , 1 , 2 , 3, are yet to be determined. 

We consider the "time-varying boundary conditions" version of the Roll 
Beam Bending Equations: viz. (4. 4), and seek a solution in the form of the 

sum 

u J(t,s) + D^(t,s) . 

It follows that we must have 

PA(U°(t,s) + D^t.s)) + EI <J u°""(t,s) = 0, 0 < s < L, 

from which we see that: 

z ( t , s) = -PAD^t.s) 

PAu£(t,s) + EI ( u°""(t,s) = -PAD^t.s) 
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Hence 

t 

u°(t,s) = / G.(t-a; s; s') (-PA D,(a,s') da 

0 9 ^ 


= (-PA) jj G^t-a; s; s') / £ a k (a) ( s ') k da ds' 


Hence the total solution is: 


3 k 3 t L 

I a k (t)s - (PA) l j J G, (t-a; s; s') a, (a) (s') k da ds’ (6.1) 

0 0 0 0 9 


The coefficient functions (a k (»)} are then to be determined from the time- 
varying boundary conditions, as we shall describe below. 


Pitch Beam Bending Equations 

Here we can write the total solution analogously as: 


3 3 t L 

l b (t)s - (PA) l j j G fl (t -a; s; s') b. (a) (s’) k da ds' (6.2) 

0 0 0 0 


where again the coefficient functions b^(») are to be determined later from 
the time-varying boundary conditions. 


Yaw Beam Torsion Equations 

Here we seek the solution in the form: 

uj(t,s) + D^(t,s) 


where 


o 



(t,s) 


satisfies: 
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P V^( t,S ) ~ GI i|; U ^"^ t » s ^ = z ( c » s > » 0 < s < L , 0 < t 

u ^' (t,0+) = u°'(t,L-) = 0 

and 

(t) s + C^(t) • 

Since we must have that: 

p y*> (t » s) + ■ G v> ,,(t ’ s) = ° » o < s < l , 

it follows that 

z(t,s) = -PI^(t,s) . 

Hence the solution is: 

t L 

C (t) + C (t)s - PI / / G.(t-a; s; s’)(C (a) + C (a)s') da ds' (6.3) 

v 0 0 v x 

where the coefficient functions Cg(t) and C^Ct) will need to be determined 
by the time-varying boundary conditions, as we shall show. 

To accommodate proof-mass conditions we will need to break up the inter- 
val [0,L] into three sub-intervals [0 , s 2 ] , [s 2 ,s 3 ] and [s 3 ,L] and choose 

different sets of the (a k (t)}, (b k (t)}, (c k (t)} coefficient functions 

as we shall illustrate later (see Section 9). 

7. FURTHER REDUCTION OF THE BOUNDARY CONDITIONS 

To proceed further with any of the methods of solution it is necessary 
to specialize the boundary functions f^ ± (t), f Q ± (t), i (t), g Q ± (t), 

j:(t), i = 1,2, 3, 4 to their given values. It is convenient to introduce 
some vector notation for this purpose. Let 
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^(t) = 


w 4 (t) = 


uj(t,0+) 

u^(t,0+) 

u^(t,0+) 

u'(t.L-) 

u^(t,L-) 

V t,L " } 


(7.1) 


(7.2) 


These are the angular velocity vectors of the shuttle and the antenna respec- 
tively. Let 


8 x (t) = 


8 <f>,l 

(t) 

s e,i 

(t) 

g >M 

(t) 

8< })»4 

(t) 

g e,4 

(t) 

g ^,4 

(t) 


g 4 (4 > " 


and let the force applied at reflector center of mass be 

,T 


Then 


F = [F ,F ,0] 
r x’ y’ J 

g 1 (t) = -(1^ +a) 1 ®I 1 w i ~ Mj^(t) - M^t)) 

g 4 (t) = -(I 4 W 4 +(o 4 ®I 4 co 4 - M 4 (t) - r®F r (t)) - m 4 r®? 4 


where 


h + m 4 


-r r 
x y 


-r r 

x y 


x 

0 


0 

0 

2 2 
r +r 
x y 


(7.3) 


(7.4) 


(7.4a) 

(7.5) 

(7.6) 


(7.6a) 
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and 


Let 


-y L -> 

u 0 (L-) 

z(L-) 


= coordinates of beam tip. 


r = 


(r ,r ,0) 

x y 


coordinates of reflector center of mass 
with respect to axes through beam tip. 


f L (t) = 


£ e ,l (t > 


(7.7) 


f 2 (t) = 


£ .K2<‘> 

f 9 > 2 (t) 


(7.8) 


f 3 (t) = 


f 6,3^ ^ 


(7.9) 


Then 


f 4 (t) = 


f e,4 (t) 


(7.10) 


f x (t) = 


m u (t,0+) 

1 <|) 

m 1 u Q (t,0+) 


(7.11) 


f 2 (t) = 


m 2 ^ ( t ’ s 2 )+ m 2%,2 

m 2 ii e (t ’ s 2 ) +m 2^e,2 


(7.12) 


f 3 (t) = 


m 3^<J)( t * s 3 ) + m 3^<j>,3 
m 3 U e (t » s 3 )+m 3 A ‘e,3 


(7.13) 


f/, (t) = -m. 


1 0 r 


x 


0 1 r 


*u e (L-) 

v l -> 


-F 


(7.14) 
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8 . METHOD II CONTINUED: ABSTRACT FORMULATION 


In this section we continue with method II, and develop an abstract 
(wave-equation) formulation of the total problem. First we begin by substi- 
tuting for the boundary functions as in the previous section. We obtain then 
the following ensemble of partial differential equations and ordinary differential 
equations . 

0 < s < s_ 


PA n^(t ,s) + El^u^'Ct, s) = 0 1 
PA u Q (t,s) + EI Q u£"(t,s) = 0 J 


S 2 < S < S 3 
s 0 < s < L 


( 8 . 1 ) 


PI ^ U ij/ t,S ^ GI ^ U ^ t,S ^ 0 » 0 < s < L 


( 8 . 2 ) 


E V5 ,(t ’ 0+) - -nyyt,o + ) 

EIqUq ' ( t , 0+) = - mi u e (t,CH-) 


(8.3) 


ElgUg" (t,L-) 


= m. 


1 0 r 


0 1 r 


V 1 ^ 

a 0 (L -> 

% (L -> 


-F y (t) 

F x (t) 


(8.4) 


E V^ (t ’ 0+) 

EI 0 Ug (t , 0+) 

G y; (t ’ o+) 


= I.M, + 


- M 1 (t) - M D (t) 


(8.5) 


EI e u J(t,L-) 

c y-J <t,I - ) 


-(I 4 OI 4 + w 4 M,(t) - roF r (t)) - ( 8 . 6 ) 
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EI <j)( u ^" (t, s 2 +) - u"' (t,s 2 -)) = -(m 2 u^(t, s 2 ) + m 2 ^(j, >2 (t)) 

EIq(uq” (t,s 2 +) - ug’(t,s 2 -)) = -(m 2 u Q (t,s 2 ) + i^Aq 2 (t)) 

| (8.7) 

EI 4>^ u ^>* ( t * s 3 + ) u ^" ( t » s 3 - )) = “( m 3 u <j > ( t » s 3^ + m 3^(|),3^ t ^ 

EI 0( U 0 M (t» s 3 +) -Uq' ( t,s 3 -)) = ~(m 3 u 0 (t,s 3 ) + m 3 ^0 >3 (t)> • J 

Let D denote the class of 3x1 functions u(*): 

V s > 

U (s ) = U 0 ( s ) 

V s) 

such that 

V u i’ u 'i> U J' e l 2 [0,L] 

and u ^' has ^-derivatives in [0,s 2 ], [s^Sg] and [s 3 ,L] ; 

u 0 » u q > u 0 » u 0 e 

and u£' has ^-derivatives in [0,s 2 ], [s 2 ,s 3 ] and [Sg,L] ; 

V u r U J € L 2 [0 ’ L] • 

Introduce the following inner products on D: 


L 



L 


L 

/ i 
0 

V s) V s) ds 

+ 

/ u Q (s) Vq ( s ) ds 

+ l > (s) V s) ds 

+ 

u 4> (0+) V 0+ > 

+ 

U 0 (0+) Vg (0+) 

+ 

y 0+) y 0+) 

-h 

u <j, (L_) v 4> (l - ) 

+ 

Uq(L-) Vq(L-) 

+ 

y ^ > (L - ) 

+ 

u* (° + ) v- (Of) 

+ 

U 0 (0+) Vg (0+) 

+ 

U 4, (L_) v i (L_) + u 0 (L_)v 0 (L_) 

+ 

u <j) (s 2 ) y s 2 } 

+ 

Ug (s 2 ) Vg(s 2 ) 

+ 

u <p (s 3 ) y s 3 } + u 0 (s 3 )v 0 (s 3 ) * 
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Complete the space under this inner product. It is readily seen that 

3 14 

the completed space denoted X will be L^COjT) *R . Define the operator 
A on D by: 

y = Ax 

where 

x = y 


v> 

E Vi" V) 


u e ( *> 

E Ve" < ' ) 


y> 

- G V *< •> 


u 4»< 0+ > 

EI^uJ' (0+) 


Ug(0+) 

EIqUq' (0+) 


U (f)( L_ ) 

- E 1 *“J’ < E -) 


Ug(L-) 

-EI e u e" 


u'(°+) 

-EI^(0 + ) 


Uq( 0 +) 

-EIgUg(0+) 


v° +) 

- c V* <0+) 


u^CL-) 

E VJ<L-> 


Uq(L-) 

EI e u e< L_ ) 


V L-) 

G Vi< L -> 


V s 2> 

EI * (u i ,,(s 2 +> 

- yy» 

u 0 (s 2 ) 

EI 0 (uJ" (s 2 +) 

“ U 0 M (s 2 - )) 

u 4 . (s 3 } 

EI 4,(“i' <s 3+) 


u 0 (s 3 > 

EI e < u 8"( s 3 + ) 

- u' ,, (s 3 -)) 


Then for u,v in D, we note that u^, v^, u 0 , v Q , u^, v^, u^, are 
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continuous in any closed subinterval of [0,L], and can integrate by parts 
to obtain: 


/ ds = u- <s 2 -)v 4 ,(s 2 ) - u-CO+jv^W) - u”(s 2 -)vj(s 2 ) 


+ u''(0+)v '(()+) + J u;;( S ) v ;;(s) ds , 


4 .'“' •<)>' 


I U™(s) v^Cs) ds 
S 2 


u i” (s 3- )v 4) (s 3 ) - u J ,(s 2 +)v <|) (s 2 ) - u J (s 3" )v i (s 3 ) 

S 3 

+ u < j ) ’(s 2 +)v ( J(s 2 ) + / uj(s) v£(s) ds , 


J U J"( S ) ds 

S 3 


u i" (L ->y L -> - - y L ->y L -> 

+ u i'< s 3 +)v i (s 3> + /V> ^ ds > 


with similar expressions for the functions Uq(»), Vg(»). Also 


L L 

/ u"(s) v (s) ds = u!(L-)v (L-) - u!(0+)v (0+) - / u'(s) v!(s) ds 


Q v ^'“ y ~~ ” J q v 4> 


From these relations it readily follows that (we omit the algebra) 

L L L 

[Au,v] = / u^(s) v^(s) ds + / Ug(s) Vg(s) ds + / u£(s) v^'(s) ds 


In particular, for u in D: 


[ Au , u ] = / (u'(s) 2 + u"(s) 2 + u"(s) 2 ) ds 


and is thus 


> 0 
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It readily follows from these calculations that A is a closed self-adjoint 
and nonnegative definite operator, and of course the domain of A is dense 
in x by the very construction of x. 

The set of equations (8.1) - (8.7) can then be expressed as an "abstract 11 
or "function space" equation, with x(t) taking values in X (see [3] 
for the general theory) 

Mx(t) + Ax(t) + Bu(t) + FN(t) + K(x(t)) = 0 (8.8) 

where M is the 17x17 matrix specified by: 


M = {m, } 

J 


terms are zero 

except 

m i,i 

= 

PA 

m 2,2 

= 

PA 

m 3,3 

= 

PI , 

m 4,4 

= 

m i 

m 5,5 

= 

m i 

m 6,6 

= 

m 4 

m 7,7 

= 

m 4 

m l3,6 

= 

m 6,13 

m 13, 7 

= 

m 7, 13 


- m, r 
4 x 


= m, r 
4 y 
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00 

00 

a 

m 8 , 9 

m 8,10 

m 9,8 

m 9,9 

m 9,10 

m 10,8 

m 10, 9 

m 10,10 

“ll.ll 

m ll,12 

m ll,13 

m 12,ll 

m 12,12 

m 12,13 

m 13,ll 

m 12,13 

m 13,13 


m 14,14 

= n 2 


m 15,15 = m 2 

m 16, 16 = m 3 

m 17, 17 = m 3 * 

We assume that M is nonsingular. Pictorially, M takes the form: 
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We note that (because it is of importance) I can be taken to be diagonal, 
the off-diagonal term being small. The "control" u(t) is 10x1: 


u(t) = 


M x (t) 

M 4 (t) 

m 2%,2 
m 2 ^ 0,2 
m 3^4> , 3 
m 3^e,3 


and B is correspondingly a 17x10 constant matrix given by 


B = 


7X10 


10X10 

(O 7 X 10 denotes 7 x 10 zero matrix) 
(IlOxio denotes 10 x 10 identity matrix) 


N(t) is the noise disturbance which is 3 x 1 , originally denoted in ( 7 . 5 ) 
by M^(t) so that F is 17 x 3 : 


F = 


7x3 


3x3 


7x3 
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Finally K(x) is a nonlinear function of x(t) given by 


K(x(t) ) 


(8.9) 


cj ® 1^ 

( 0 . ® 1 , 0 ), 

4 4 4 


The most important thing to note about the function K(x(t)) is that 


Indeed 


[K(x(t) ) , x(t) ] = 0 


[o^®!^, 0^] = 0 


W 4 ] = 0 


The abstract formulation has many advantages beyond the immediate one of 
providing a succinct statement of the essence of the problem. First 
the total energy in the system is 


E(t) = j [Mx(t), x(t)] + \ fAx(t), x(t)] , (8.13) 
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the time-derivative of which is: 

E(t) = [Mx(t) + Ax(t), x(t) ] 

which, using (8.8) 

= - [Bu(t) , x(t)] - [K(x(t) ) , x(t) ] - [FN(t), x(t)] 

= - [Bu(t) , x(t) ] - (FN(t) , x ( t) ] 

since 

[K(x(t)) , x(t)] = 0 . 


Hence, if we ignore the noise term for the moment, we see immediately that 
the system can be "stabilized" in the sense that the system energy E(t) 
decreases (or does not increase) by taking the feedback control: 


M 1 (t) 


^(t) 

M 4 (t) 


<o 4 (t) 

*• 

= P 

y^y 

/— s 
■U 

Sm/ 

CM 

CD 

<3 

CM 

£ 


u 0 (t,s 2 ) 

m 3V3 (t) 



m 3^0 ,3^^ 


u Q (t, s 3 ) 


where P is positive definite. In particular it is enough to take P 


the form 


P O 

6x6 u 6x4 


°4x6 °4*4 


where the 6x6 matrix P 6xg is positive definite. We can thus dispense 
with the "proof-mass" controllers and use only "feeding back" appropriately 


the angular velocities. 
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To proceed further with solving the abstract equation (8.8), let 


Y(t) = 


x(t) 

x(t) 


Then we can rewrite (8.8) in the usual "state space" form as 


Y (t) = AY(t) + flU(t) + K( Y(t)) + FN(t) 


where 


A = 


0 I 

-M _1 A 0 


Let 


Then 


5U = 


Y = 


H Y) = 


FN(t) = 


0 

-BU 


0 

-K(y 2 ) 

0 

-FN(t) 


We next introduce the energy inner product (see [3] for this) 


[Y,Z] e = [My 2 ,z 2 ] + [v^y^ , ] 


where 


Y = 


Z = 


(8.16) 


(8.17) 


(8.18) 


(8.19) 


( 8 . 20 ) 
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Then in this inner product 



so that A is dissipative: 

[AY, Y] e + [Y, 4Y] e < 0 . 

Moreover A has a compact resolvent and the spectrum of A consists entirely 
of the "point" spectrum (eigenvalues) . Now the eigenvectors of A given by 

AY = XY (8.21) 

or, using (8.18), we have 

y 2 = Xy l 

-M 1 Ay 1 = Xy 2 = X 2 y x 

or, 

Ay 1 = -X 2 My 1 (8.22) 

Since A is self-adjoint and nonnegative definite and so is M it follows 
that we can express the eigenvalues as 
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n 


n 


+i y en 
n 


n 

-i u e 
n n 


where 


Ae « y Me 
n n n 


and the (e } is an M-orthogonal basis in X. In other words we are 


n 


guaranteed to find functions e^Cs), 0 < s < L such that, writing 





* n < s ) 


e = 


e (s) 


n 


n 




V*> 

<f£"(s) 

2 2 
= Ya y 
9 n 

V s) ] 





CD 

3 Z 

CO 

2 2 
= Vn 

9 „< s > J 

~K (s) 

2 2 
= Vn 

V s) 


in the open intervals 
(Oj s 2^ * ^ s 2 ,s 3^ * (SjjL) 


0 < s < L 


El . <f) ,M (0+) 

<p n 


EI ft O'" (0+) = 

u n 


El <J>'" (L-) 
9 n 


EI e 9™ (L-) 


m, V * (0+) 

Inn 


m. p 0 (0+) 

I n n 


= -y 


'^(Vn^ + m 4 r x^n (L-) 
m 4 0 n (L_) + m 4 r y V L_) ) 


(8.23) 


(8.24) 


(8.25) 
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EI «, ♦”(0+) 


<j>;(°+) 

EI 0 0"<O+) 

" T 1 

0;(o+) 

GI * K (0+) 


^ n (o+) 


l 


EI c» 




+ ;a-) 


0 

EI 0 

9 n< L -) 

= V 

A 

*4 

%< L ’> 

+ 

0 


^( L -) 


\ 

*n (L -> 


m.r 4> (L-)+m,r 0 (L-) 
4 x n 4 y n 


E y*n (s 2 +) - C (s 2-» 
EI 0 (0 n (s 2 +) - 6 n < s 2-» 

E V*n (s 3 +) - ♦» (s 3- }) 

El (0'” (s +) - 0'" (s -)) 

0 n o no 

Moreoever for distinct eigenvalues 

[Ae , e ] = 0 = 

n m 


= m 2 •‘n W 


= m y 0 (s ) 
2 n n 2 


= m y i(sj 

o n n o 


= m y 0 (s ) 
o n n o 


n 


y 


m 


[Me , el 
n m 


and we may orthonormalize so that [Me ,e ] = 1. Note that we do not need 

n n 

to verify whether any matrices are nonsingular as in [2] for exmaple! We shall 
go into the details of actually calculating the eigen functions in Part 2. 

Next we can write a "modal expansion" for the solution of (8.8). Thus let 

CO 

= l a (t) e , (8.26) 

0 

2 

where we take the eigenvalues y^ in increasing order, the smallest eigenvalue 
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being y Q . Substituting (8.26) into (8.8). we obtain 

00 CO 

1 (5 n (t) Me n + y n 5 n (t) Me n ) = l h n (t) e n <8.27) 

where we take 

h (t) = [h(t) , Me ] 

n n 

h(t) = -FN(t) - Bu (t) - K(x (t) ) . 

Hence 

[ 5 n (t) + V n (t)] = h n (t) [e n ,e n ] 

so that 

t y n (t-a) 

a n (t) = J sin h n (a) da [e ,e ] (8.28) 

0 P n 

taking initial conditions to be zero. Hence we can write (in terms of the 
operator valued Green's Function): 

t 

x(t) = J G(t-a) h(a) da (8 29) 

0 

where G(t) is defined by 

00 sin y t 

G (Ox = l - — — [e ,e ] [x,Me ] e 

0 M n n n n n 

for each x in X. To allow for the nonlinearity we must next solve the 
integral equation: 
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t t t 

x(t) + J G(t-a) K(x(a)) da + J G(t-a) FN(a) da + / G(t-a) Bu(a) da 
0 0 0 

= 0 

We shall examine this equation in more detail in Part 2. 

9. AN EXAMPLE OF AN EXPLICIT SOLUTION OBTAINED BY THE BOUNDARY INPUT METHOD 

In this section we develop one explicit solution to the Beam Equations 
by means of the Boundary Input technique outlined in Section 6. 

The simplest such solution is obtained by setting, using the notation 
of that section, 

fi^Ct.s) = 0 ' 

ij Q (t,s) = 0 > (9.1) 

V t,s) = 0 j 

so that the total solution (6.1), (6.2), (6.3) can be expressed very simply as 

3 k 

u <j) (t > s ) ■ l a k (t) s 

3 k 

u 0 (t,s) = l b (t) s 

yt, s ) = c 1 (t)s + c 0 (t> 

All we need to do now is to determine the coefficients in (9.2) (taking into 
account (9.1)) by substituting into the equations (8.3) through (8.7). Since 
we are only interested in exhibiting one particular solution, we shall simplify 
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the problem by neglecting the proof-mass conditions. We shall also neglect 
noise, so that 

MpCt) = 0 . 

Next we demand a solution in which the angular accelerations are zero: 


(S 1 = 0 


% = o 


Condition (8.3) yields that 


a 3 (t) = 0 


b^ (t) = 0 . 


Condition (8.4) is satisfied by taking 


y t ) ■ 

/-N 

rt 

v—' 

= 0 

• 

first that 

yt) 


w 4>,i 

“i = 

yt) 

= 

w e,i 


c Q (t) 


CO. . 
^,1 


(9.3) 


Since (0^ is constant, this implies that 


a l (t) = tW (f),l + a i 

b 1 (t) = tw e } i + b i 

c 0 (t) = tco^ + c Q 


(9.4) 
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W <M = 2a 2 (t)L + a-^t) 


Jq 4 = 2b 2 (t)L + b-^t) 


%4 = d l (t)L + d 0 (t) 


a 2 (t) = 


/. ^A. 1 ) 


b 2 (t) = 


( “0.4 - Vl> 


(w, . - 0) , , ) 
CjCt) - »' 4 - I 1 ' 1 


Setting the initial "position 11 conditions: 


u cj)(0» s) 

Uq (0,s) 
u^(0,s) 


^4>(o» s) = 
D 0 (O,s) = 

D^(0,s) = 
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So that, in particular, in (9.4) 


a l b l c 0 ~ 0 * 


a 2 (0) = b 2 (0) = Cl (0) = 0 , 


thus determining completely the functions a 2 (t), t> 2 (t), a^t), b^(t), 

c^(t), c Q (t) . To satisfy (8.5) and (8.6), we define the control moments 
M^(t) and M^(t) by 


t (*v 


( 03 . . - 03 , - ) 

9,4 _ . <5,1 


(“ fl , - oj ) 

M x (t) = t(EI 0 ) - - U > 1 - - ( o 1 ® 1^ 


t < G V 


(03 . . - 03 , , ) 

<5,4 T|3,l y 

L 


M 4 (t) = -M^t) 

Thus the final solution is: 



V t,s) 

= t 

(< V 

s 2 

2L S 

+ 

tW 4>,l 


Ug(t,s) 

- t 

( “0,4 

~ “M* s 2 

2L S 

+ 

tw e,i 

♦ 

V t,s) 

= t 

<V 

'Pill s 

L s 

+ 

tc V 


where a^, b Q are arbitrary constants, and can be chosen to satisfy initial 
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time-derivative conditions for u^(t,s) and Ug(t,s). 
energy in the system, defined by (8.13) is a constant, 
are 


u^(0,s) 

Uq(0,s) 

u^(0,s) 


(M *.< ' "♦.i ) j + 

2L S + 

(6) e,4 ~ VP 2 

OT ® * 


(&>, 




-CO. 


+ 


Note that the total 
Also, the initial "rates" 

w 4> ,1 s + a 4» 

w e,i s + a e 


J 


which are of course nonzero. The control is chosen, roughly speaking, to 
maintain this rate. Our aim here has been limited to obtaining one explicit 
solution. In the forthcoming work we shall enlarge further on this technique. 


CONCLUDING REMARKS 


Several techniques of solution for the partial differential equation 
formulation of the SCOLE control problem have been described which should 
facilitate conceptual understanding as well as numerical computation of feed- 
back laws. Additional work is required before the relative merits of the 
various approaches can be evaluated. 
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